Knotting of linear DNA in nano-slits and nano-channels: a numerical study 
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The amount and type of self-entanglement of DNA filaments is signficantly affected by spatial confinement, 
which is ubiquitous in biological systems. Motivated by recent advancements in single DNA molecule exper- 
iments based on nanofluidic devices, and by the introduction of algorithms capable of detecting knots in open 
chains, we investigate numerically the entanglement of linear, open DNA chains confined inside nano-slits. The 
results regard the abundance, type and length of occurring knots and are compared with recent findings for 
DNA inside nano-channels. In both cases, the width of the confining region, D, spans the 30nm- l/^m range 
and the confined DNA chains are 1 to 4/im long. It is found that the knotting probability is maximum for slit 
widths in the 70-IOOnm range. However, over the considered DNA contour lengths, the maximum incidence of 
knots remains below 20%, while for channel confinement it tops 50%. Further differences of the entanglement 
are seen for the average contour length of the knotted region which drops significantly below D ~ lOOnm for 
channel-confinement, while it stays approximately constant for slit-like confinement. These properties ought 
to reverberate in different kinetic properties of linear DNA depending on confinement and could be detectable 
experimentally or exploitable in nano-technological applications. 



I. INTRODUCTION 

Like other forms of entanglement, the abundance and type 
of knots in equihbrated DNA molecules depends on both in- 
trinsic and extrinsic properties. The former include the chain 
contour length, the bending rigidity and chirality while the lat- 
ter include spatial constraints such as confinement in narrow 
spaces 1(131 [T7ll22l. 

A notable example of spatially-confined DNA is offered 
by viral genome packaged inside capsids[13|. Typically, the 
viral DNA contour length exceeds by orders of magnitude 
the capsid diameter, resulting in a near crystalline density of 
the packaged genome |7 9, 10, 23 1. Over the years, sev- 
eral numerical studies have accordingly tried to understand 
not only how the DNA filament is packaged inside the virus 
C] El 13 H] [B] |23l 1211 but especially how it can be ejected 
into the host cell through the narrow capsid exit channel with- 
out being jammed by self-entanglement |2] [141 US 123 ■ A 
solution to this conundrum was proposed in ref lfT4ll which 
reported that the ordering effect of DNA cholesteric self- 
interaction [4, 29J is responsible for keeping the entanglement 
at a minimum and compatible with an effective ejection pro- 
cess. 

The effects of spatial constraints on DNA self-entanglement 
and the possible implications on DNA condensation, pack- 
aging and translocation have been systematically addressed 
only recently, largely because of the introduction of suitable 
nano-devices and micro-manipulation techniques that allow 
for probing the properties of few confined molecules at a 
time |l3l|5l|25]|26,i0,i2J. 
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In such contexts a still largely-unexplored research avenue 
is the characterization of the occuiTence of knots in open, 
linear DNA molecules. In fact, theoretical and experimen- 
tal studies of knot occurrence have largely focused on equili- 
brated chains where knots are trapped by a circularization re- 
action which ligates the two chain ends, thus forming a ring. 
The topology of such rings, is clearly maintained until they 
open up, and therefore their knottedness is well-defined. 

This is not the case for open chains, where non-trivial en- 
tanglement cannot be permanently trapped because of the 
two free ends. Yet, we are all familiar with the fact that 
knots in open chains can be long-lived and can affect vari- 
ous physical and dynamical properties of polymers. In par- 
ticular in lab-on-chip experiments the presence of knots in 
linear DNAs may interfere with the confinement elongation 
process of the molecules, an essential step for the detection 
of protein-DNA interactions |36| and also towards genome 
sequencing by pore-translocation If27l 1371 . These considera- 
tions have stimulated a number of efforts aimed at suitably ex- 
tending the algorithmic notion of knottedness to linear, open 
chains 121] [33l |35l . 

Building on these theoretical advancements and motivated 
by the upsurge of DNA nano-manipulation experiments IS] |5l 
|251 l26l [30l [32l here we report on a numerical study of the 
knotting properties of linear DNA chains confined in nano- 
slits and nano-channels. 

The investigation is based on a coarse-grained model of 
DNA and is a foUowup of two recent studies of the metric 
and entanglement properties that we carried out for closed and 
open chains in nano-slits and nano-channels f T9ll20l . Specifi- 
cally, the properties of knotted open chains confined in slits is 
reported here for the first time and is compared with the earlier 
results for channel confinement. 
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II. METHODS 

The model 

Hereafter we provide a brief, yet self-contained, description 
of the coarse-grained DNA model and simulation and numer- 
ical techniques used to characterize the topological properties 
in slit- and channel-like confining geometries. 

By following the approach of ref. |28|, linear filaments of 
dsDNA are modeled as semi-flexible chains of N identical 
cylinders. 

The cylinders diameter is d — 2.5nm, corresponding to the 
dsDNA hydration diameter, and their axis length is equal to 
b = lOnm, i.e. a fraction of the nominal dsDNA persistence 
length, Ip ~ 50nm. 

A chain configuration is fully specified by the location in 
space of the cylinders axes, ti, t2,tN and, in the uncon- 
strained case, its energy is given by the sum of two terms; 



We considered chains of N ==120, 240, 360, 400 and 480 
cylinders, corresponding to contour lengths Lc = Nb in the 
range 1^.8/j,m. Across the various values of fj, we collected 
^ 10^ uncorrected configurations of chains that could be 
accommodated inside channels or slits with width D in the 
40 nm-l/im range. Notice that, because of excluded vol- 
ume effects, the minimal width achievable by slits and chan- 
nels in presence of a linear chain is d. It is therefore conve- 
nient to profile all properties in terms of the effective width 
D,s^D- d. 

Chain size and knot detection To characterize the average 
size of the chain we consider its root mean square radius of 
gyration 
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The first term accounts for the excluded volume interaction of 
the cylinders. It is equal to "infinity" if two non-consecutive 
cylinders overlap and zero otherwise. The second term is the 
bending energy, Et = -KgT'f ^ ' ^+i with T = 3G0K 
being the system temperature and Kb the Boltzmann con- 
stant. It is assumed that the DNA is in a solution with high 
concentration of monovalent counterions, so that its electro- 
static self -repulsion is effectively screened. Because non-local 
self-contacts of the DNA molecule are infrequent for two- and 
one-dimensional confinement we neglect both desolvation ef- 
fects and the DNA cholesteric interaction 11141 [311 . Finally, 
because linear chains can effectively relax torsion, the DNA 
torsional rigidity is neglected too. 

The confinement of the DNA inside slits is enforced by re- 
quiring that the chain maximum span perpendicular to the slit 
plane , A, is lower than a preassigned value, D. Likewise, for 
channel confinement, it is required that the maximum calliper 
(diameter). A, measured perpendicularly to the channel axis 
is smaller than D. 

The conformational space of confined chains was explored 
by means of a Monte Carlo scheme employing standard lo- 
cal and global moves (crankshaft and pivot moves). Fol- 
lowing the Metropolis criterion, a newly-generated trial con- 
formation is accepted or rejected with probability given by 
min(l, exp[— /iA)/iirBT]. In the latter expression rep- 
resents an auxiliar parameter that couples to the chain span (or 
calliper size), A. Accordingly, by using different values of /i it 
is possible to bias the sampling of the configurations towards 
configurations with different average values of A. Next, be- 
cause the biasing weight is set a priori, it is possible to re- 
move it by using thermodynamic reweighting techniques and 
recover the canonical expectation values of the observables 
of interest. Advanced sampling and reweighting techniques 
(which are reviewed in detail in ref. f\n\) is adopted here be- 
cause a direct enforcement of the geometrical constraints in 
the Monte Carlo sampling would be inefficient due to the high 
Metropolis rejection rate. 



where a runs over the three Cartesian components of the posi- 
tion vector Yi of the center of the z-esim cylinder of the chain, 
f = X^iLi is '^he center of mass of the chain and the 
brackets (•••) denote the canonical average over chains in 
equilibrium an confined in slits or channels. 
ProflUng the knot spectrum 

The entanglement properties of the confined model DNA 
filaments were characterised by establishing the knotted state 
of the open chains and by measuring the knot contour length. 

From a mathematical point of view only circular chains 
have a well-defined topological knotted state since it cannot 
be altered by distorting or changing the chain geometry as 
long as the chain connectivity is preserved. To extend the con- 
cept of knottedness to open chains, it is therefore necessary to 
close them in a ring with an auxiliar arc 1211 [33l . The knot- 
ted state of the closed ring is then assigned to the open chain. 
The auxiliary arc must clearly be suitably defined to ensure 
the robustness of the topological assignment; in particular it 
must avoid interfering with the self-entanglement of the open 
chain. To this purpose we adopted the minimally-interfering 
closure scheme introduced in ref. |33|. 

The position of the knot along the chain is next established 
by identifying the shortest chain portion that, upon closure 
has the same knotted state of the whole chain. To minimize 
the chance of detecting sUpknots ETl it is also required that 
the complementary arc on the closed chain is unknotted. 

Fig.[T]illustrates two knotted configurations of 2.4/im-long 
open chains confined inside a slit and a channel. The knotted 
portion of the chains is highlighted. 



III. RESULTS AND DISCUSSION 

The metric properties of linear DNA chains for various de- 
grees of slit-like confinement were systematically addressed 
in ref |20|. Such study indicated that for increasing confine- 
ment, the two principal axes of inertia of linear molecules first 
orient in the slit plane and next grow progressively as the chain 
spreads out in a quasi-two-dimensional geometry. 
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FIG. 1. Linear knotted chain of length Lc ~ confined inside a slit and a channel. In both cases the confinement width is D 

and the chains accommodate a trefoil knot (highlighted in red). 



70nm 



The interplay of the increase of the chain size projected in 
the slit plane, R\ | , and the concomitant decrease of the trans- 
verse size, R±, results in the non-monotonic behaviour of Rg, 
as illustrated in Fig. |2] The data shows that Rg attains a mini- 
mum at an effective channel width, D*, that is slightly larger 
then the average radius of gyration of the unconstrained, bulk 
case (marked by the black dashed line). 

For comparative purposes, in the same figure it is shown the 
average radius of gyration of equally-long linear DNA chains 
confined in channels. In this case too one observes the non- 
monotonic dependence of Rg on the width of the confining 
region, which attains its minimum at a channel width, D* 
slightly larger that the slit case (red dashed line). However, 
the increase of Rg past the minimum is much more dramatic 
than for the slit case. This reflects the fact that the chain can 
only elongate in only one direction rather than in a plane. 
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FIG. 2. Average radius of gyration, Rg, of linear DNA with contour 
length Lc = 4.0/im inside a slit and a cylindrical channel of effective 
width -Deff. Black and red dashed lines mark the location of the 
minimum of Rg respectively for slits and channels. 

A relevant question regards the extent to which the di- 
mensionality of the confining region (slits or channels) and 
its width can affect the incidence, complexity and size of 
knots in linear chains. We accordingly applied the minimally- 
interfering closure scheme to establish the knotted state of 
equilibrated linear chains and to locate the knot along their 
contour. 

We first discuss the overall incidence of non-trivial knot 
topologies. The results for slit-confined chains are shown in 
Fig.[3]^a). As expected, for each fixed value of D^ff the knot- 
ting probability depends strongly on Lc lil7J . In fact, going 



from Lc = 1/im to 4.8/im the knotting probability increases 
by one order of magnitude. By comparison, the knotting prob- 
ability variations on Doff (at fixed Lc) are smaller, though no- 
ticeable. More importantly, the knotting probability displays, 
as DcH decreases, a non-monotonic behaviour with a maxi- 
mum enhancement peak at a width, Dc that falls within the 
50-100 nm range. 

In particular, the knotting probability varies by a factor of 2 
going from the unconstrained case D ^ Ifim to Dc ^ 80nm. 
As consistently indicated by analogous results of slit-confined 
rings 1 20 1, this knotting enhancement should be measurable 
experimentally by circularizing dsDNA molecules with com- 
plementary single-stranded ends inside slits. 

As shown in Fig. [3|b), the confinement-induced enhance- 
ment of non-trivial knots is even stronger for the channels 
case: at the largest contour length, Lc = 4.8/im, the prob- 
ability peak value is about 10 times larger than the bulk one. 
One may anticipate that the maximum knot enhancement is 
attained for the same channel width, D*, for which Rg is min- 
imum ( i.e. at the highest value of the overall chain density). 
This is, however, not the case since Dc is about one third of 
_D* at all considered chain lengths. 

Besides the overall knotting probability, it is of interest to 
examine in which proportion the various types of knots con- 
tribute to the overall knot population. The results for 4.8/im- 
long chains inside slits are shown in Fig. [4] One notices that, 
at all explored values of Dcs, prime and composite knots with 
up to six crossings in their minimal diagrammatic representa- 
tion account for at least 90% of the observed knot population. 
In particular the simplest knot type, the 3i or trefoil knot, is 
by far the most abundant. For channel confinement the pre- 
dominance of simple knots is even stronger In particular the 
peak probability of the trefoil knot is, at _Doff = Dc about 
23%, which is four times larger than for the unconstrained, 
bulk case. 

The result is noteworthy for several reasons. Firstly, at 
variance with the case of three-dimensional isotropic con- 
finement (cavity, capsids) of DNA rings ifTSll the knot spec- 
trum of chains in slits and channels is dominated by the sim- 
plest knot types. This fact was recently established for closed 
chains in slits and channels and for open chains in channels 
only |J9u 20|. The present results for open chains inside slits 
therefore complete the overall picture in a consistent way. 
Secondly, the percentage of simple knots (and in particular 
of the simplest one, the trefoil) found in two-dimensional con- 
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FIG. 3. Percentage of knotted linear DNA chains inside slits (a) and cylindrical channels (b) of transversal effective size 7?cff- Different curves 
refer to different DNA contour lengths. The estimated relative errors on the knotting probability are smaller than 1% for Deft > 200 nm and 
about of 3 — 4% for Deft < 200 nm. The knotting data for linear chains inside channels are based on the study of ref. 1191 . 



finement (slits) is at least doubled in one-dimensional confine- 
ment (i.e. channels). 

To profile the finer characteristics of the chain self- 
entanglement, we finally analysed the average length of the 
chain portion that is spanned by the highly abundant trefoil 
knots, ^3j. The results shown in Fig. [5] show that, at fixed 
Lc, is non-monotonic on Doff for both slits and channels. 
However, one major difference is that is that the knot length 
decreases dramatically after the peak for channels while for 
slits it appears to approach a limiting value not dissimilar from 
the bulk one. 

For both types of confinement, depends strongly on L^- 
For instance, the confinement width associated with the max- 
imum knot length varies noticeably with Lc- In addition, the 
relative difference of the peak value of respect to the bulk 
case increases with Lc too (it ranges from 10% to 22% for 
slits and from 16% to 36% for the channels). 

It is interesting to examine the observed dependence of 
^3^ on Lc in connection with earlier scaling studies on 
closed rings in various conditions (unconstrained, collapsed, 
stretched etc. |^[TTl[T2][34l)- In particular, we recall that knots 
in unconstrained rings are known to be weakly localized in 
that their contour length scales sublinearly with Lc- Specif- 
ically, for trefoils it was shown that for asymptotically-long 
rings, ?3^ - L* with t w 0.65 U2J. 

Motivated by these earlier findings we analysed the data in 
Fig.|5]and rescaled them so to collapse the ^3^ curves for very 
weak confinement (i^off > SOOnm). The optimal rescaling 
was obtained for the exponent, t = 0.62 ± 0.05 for both slits 
and channels. For stronger confinement, the rescaled curves 
depart from each other By contrast with the slit case, the 
channel data show a systematic upward trend for increasing 
Lc- As a matter of fact, the portion of the curve for D^s < Df, 
shows a good collapse for i = 0.74±0.01. 

This suggests that, as confinement increases, the knotted 
subregion of the chain remains weakly localized but it further 
swells along the unconstrained dimensions. 



IV. CONCLUSIONS 



In polymer physics there is an ongoing effort to understand 
the extent to which spatial constraints affect the probability 
of occurrence, the complexity and size of topological defects 
in linear polymers. For DNA this problem has various im- 
plications both for the understanding of some biological el- 
ementary processes (such as translocation and viral ejection) 
and for the development of efficient setups for DNA nano- 
manipulation protocols (such as sorting or sequencing). 

Here we reported on a numerical study, based on ad- 
vanced Monte Carlo simulations, thermodynamic reweight- 
ing and scaling analysis of the equilibrium topological prop- 
erties of a coarse grained model of linear DNA confined in 
rectangular slits and cylindrical nanochannels. The investi- 
gation was carried out for linear chains with contour lengths 
ranging between 1.2 and 4.8/Lim and confined within geome- 
tries whose transversal dimension DcS span continuously the 
30 — lOOOnm range. 

We found that, both for slits and channels, the knotting 
probability is a non monotonic function of D^s with a peak 
that occurs at a length-dependent confinement width Dc- 
Most importantly and unlike DNA in capsids (i.e. under full 
confinement) the enhancement of the topological entangle- 
ment in slits and channels is not followed by a corresponding 
enhancement of the entanglement complexity. Indeed, despite 
that the peak knotting probability exceeds by several times the 
one in the bulk, most of the knots observed belong to very sim- 
ple knot types. This effect is particularly evident for channel 
confinement. This suggests that nano-fluidic devices based on 
this or similar one-dimensional geometries may be very effec- 
tive for producing a good population of linear DNA molecules 
with a simple knot tied in. 

Finally, by using a robust algorithm for locating knots in 
open chains |33|, we show that the typical contour length of 
the knotted region displays a non monotonic behaviour similar 
to the one observed for the knotting probability. Moreover, by 
looking at its scaling behaviour as a function of the chain con- 
tour length, it is found that for the whole range of confinement 
and both for slits and channels, knots are weakly localized. 
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FIG. 4. Knot spectrum as a function of D^h for 4.8^m-long linear DNA chains confined within slits (a) and channels (b). The knotting data 
for linear chains inside channels are based on the study of ref. 1191 . 
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FIG. 5. Average contour length of knotted arcs in linear DNA in slits (a) and channels (b) of effective transverse dimension Deff- In (c) and (b) 
the data of (a) and (b) are rescaled by L* for Ig-^ and by L" with v = 0.585 for Deff- The main panels refer respectively to t = 0.62 ± 0.05 
while for the plot in the inset t = 0.74 ± 0.01. 
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